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Abstract 

We present the results of our studies of the entanglement entropy of a superconducting system 
described holographically as a fully back-reacted gravity system, with a stable ground state. 
We use the holographic prescription for the entanglement entropy We uncover the behavior 
of the entropy across the superconducting phase transition, showing the reorganization of the 
degrees of freedom of the system. We exhibit the behaviour of the entanglement entropy from 
the superconducting transition all the way down to the ground state at T = 0. In some cases, 
we also observe a novel transition in the entanglement entropy at intermediate temperatures, 
resulting from the detection of an additional length scale. 
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1 Introduction 



Given a quantum system, the entanglement entropy of a subsystem A and its complement B is 
denned as follows: 

Sa = -Tr.4 (P^lnp^i) , (1) 

where pj± is the reduced density matrix of A given by tracing over the degrees of freedom of £>, 
PA = Trg(/?), where p is the density matrix of the system. 

The entanglement entropy is understood as an important probe of physics in various do- 
mains, and for systems at strong coupling it is looked upon as a robust tool for keeping track 
of the degrees of freedom when other traditional probes (such as an order parameter) might not 
be available. However, it is often difficult to compute the entanglement entropy in such systems, 
especially outside 1+1 dimensions. 

Two developments in the field have made the work in this paper possible. The first is that 
the entanglement entropy has a natural geometrical definition [TJ [2] (proposed but only partially 
proveij^J) in the context of gauge/gravity duals, where a wide class of strongly coupled theories in 
d dimensions can be defined holographically as dual to a theory of gravity (plus other degrees of 
freedom) in d + 1 dimensions, which is in turn ultimately embedded into a 10 dimensional super- 
string theory or an 11 dimensional supergravity background as a means of ensuring full quantum 
consistency. Of course, holography is also not fully proven, but there is a large body of evidence 
for it in numerous examples, starting with the AdS/CFT correspondence and deformations and 
generalizations thereof [111 [T2~l [13| I14j . In this paper we shall assume that holography is robust, 
and that the holographic formula (reviewed below) for the entanglement entropy is also correct. 

The second development is that some of the strongly coupled physics of interest, supercon- 
ducting phases that share various features with certain exotic phases of experimentally studied 
strongly coupled quantum systems, can not only be modelled holographically as an effective model 
of gravity plus a scalar (see e.g. ref. [IS])) not only be embedded consistently into the parent super- 
string theory and/or 11D supergravity to get access to the back reactions of the scalar dynamics 
on the geometry (see e.g. ref. [E]), but can be embedded in a manner that appears to be highly 
stafr/^] (i.e., refs. [H] and [20] have shown that the ground states of an infinite subset of the family 
of superconductors defined by the embeddings in ref. |17j are in fact unstable in maximal N = 8 

'See e.g., refs. H O H El H E3] • 

2 To be precise, the physics breaks a global, not local, symmetry, but it is close to being gauged, in a sense [15] . 
So rather than using the term superfluidity, we will continue with the common usage. 

3 Strictly speaking, the stability is typically studied in various truncations of 11D supergravity to lower dimensions. 
We mean here full perturbative stability of the ground state in maximal N = 8 supergravity in 4D. There remains the 
possibility of instabilities arising upon uplift to the full 11D supergravity, non-perturbative instabilities, and parts of 
the phase diagram being modified at higher temperatures by instabilities of the sort discussed in e.g. ref. [18] . 
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supergravity in 4D and it remains to be shown whether any of the others in the family are stable). 
In other words, there is a holographic superconductor model (presented recently in ref. [21]) that 
has a ground state that, thought of as a holographic flow [22] is a stable^] non-supersymmetric vac- 
uum of the full theory. This suggests that the complete theory of gravity plus all the attendant 
fields is without pathological physics that might obscure the lessons to be learned from it about 
strongly coupled phenomena pertinent to the superconductivity. 

These two developments come together nicely since to employ the holographic definition of 
the entanglement entropy in a study of superconductivity, we need the complete (back-reacted- 
upon and stable) geometry to perform the computation. In this paper we carry out the study of 
the entanglement entropy in holographic superconductivity for the first time using these methods, 
and uncover some very interesting phenomena. 

The entropy is holographically computed as follows [H [2]. In an asymptotically Anti-de 
Sitter (AdS) geometry, consider a slice at constant AdS radial coordinate z = a. Recall that this 
defines the dual field theory (with one dimension fewer) as essentially residing on that slice in the 
presence of a UV cutoff set by the position of the slice. Sending the slice to the AdS boundary at 
infinity removes the cutoff (see ref. |23j for a review). On our z = a slice, consider a region A. Now 
find the minimal-area surface 7^4 bounded by the perimeter of A and that extends into the bulk of 
the geometry. (Figure [T] shows examples of the arrangement we will consider in this paper.) Then 
the entanglement entropy of region A with B is given by: 

Area(7^) 

where Gn is Newton's constant in the dual gravity theory. 

t 




Figure 1: Diagram of the strip shape we will consider for region A. This is the case of a dual geometry 
that is asymptotically AdS4, and here, z denotes the radial direction in AdS4. The quantity £ sets the size 
of region A, and L is a regulator that is understood to be taken to infinity. 



In the next section, we will review the four dimensional model of gravity plus scalars and a 
4 Subject to the caveats in the previous footnote. 
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gauge field that was presented in ref. |21j . and review and re-derive the properties of the solutions 
we need. 



2 Gravity Background 

The Lagrangian of ref. [21] arises as an SO(3) x 50(3) invariant truncation of four-dimensional 
M = 8 gauged supergravity: 

£ 1 1^ ^,,„ n sinh 2 (2A) /„ o J \ / „ q 



where the potential V is given by: 



(3) 



P - -5 2 ( 6 cosh 4 A - 8 cosh 2 A sinh 2 A + - sinh 4 A ) . ( I i 



We use a different notation from ref. [2T]. We first reintroduced the dimensionful constant 87^4 
and then made the following field redefinitions: 

Note that the gauge field and the (real) scalar fields A and ip are dimensionless in this framework. 
The metric and other fields of interest are parameterized as follows: 

ds 2 = -^f( z )e~^dt 2 + % {dx\ + dxl) + % ^ , At = *(z), A = X(z) , (6) 
and the scalar ip will be set to zero. Defining a useful dimensionless coordinate: 

z = Rz , (7) 

the equations of motion can be reduced to: 

/-.W- ^f'^ ,, (8) 

„ fV\ . sinh 2 (2A)# . , 

* + (fj V = , (10) 

.,,,( X' /' 2\ i? 2 e^sinM4A^ _ 

A A "4^/dA + 16^ " ( } 
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The ultraviolet (UV) asymptotic behavior (near the AdS boundary z = 0): 



X(z) = X 1 z + X 2 z 2 + ... , 

X(z) = Xo + Aq2 2 + . . . , 

f(z) = 1 + X%z 2 + f 3 z 3 + . . . , 

*(z) = ¥ + *l* + ..- ■ (12) 

Generically we will be at finite temperature, to which there will be associated an event horizon in 
the geometry. We assume the event horizon occurs at z = zjj, and near there the fields have an 
expansion: 



A(jg) = A(°) + A«(l-~)+... 



,(o) , v (i) i - 



= U-^J + ... • (13) 

There are only three independent parameters here, which we choose to be A (0) , x {0) , There 
are three scaling symmetries of the equations of motion given by |21j : 

t->7f x t, x^X-21n 7 i, tf-^tf, (14) 
i — > 7 2 T 1 i , z — ► 7 5 r 1 z , -R ->■ 7 2 ~ 1 -R , (15) 
^->7-V, /(*)->/(«), *(^)^7*(^) J A(z)->A(z), x(«)->x(2). (16) 

Using these scaling symmetries, we can choose any value for the position of the event horizon and 
the asymptotic value of the field x( z )- We choose to fix A^ ) and x^°\ an d tune VljW to either have Ai 
or A2 be zero. These asymptotic values of the field A in the UV define the vacuum expectation value 
(vev) of charged operators in the theory that are either of dimension 1 or 2, and we correspondingly 
call them 0\ and 02- We will explicitly identify the correctly normalized relationship below. The 
UV asymptotics of the electric gauge field component *&(z) defines a chemical potential and charge 
density that will be explicitly identified below. 

Generically, a solution will have xo non-zero. To recover pure AdS, we use the scaling 



symmetry in equation (14) to shift x( z ) as: 

X(z) = x(z) - Xo , (17) 
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which can be accomplished by rescaling the time coordinate: 

t = e~ xo/2 t, (18) 

which in turn means: 

A- t = e xo/2 A t . (19) 
The temperature of the system is then given by [21] : 

-(x (0) -Xo)/2 



T = 4 J^g g 3 ~ (ei + 36cosh ( 2 A<°>) - cosh (iX^) - Sz^ (V 1 ))*) 



and the chemical potential p and charge density p go as: 



and the vevs of the two operators are defined as: 



Using the scaling symmetry in equation (16), the relevant quantities for us scale as: 



Therefore, we will use the following dimensionless quantities to examine the physics: 

JL 9±_ 9i 

y/P ' y/P ' P 

2.1 High Temperature Phase 



(20) 



e xo/2 e x °/ 2 

» = 7muf^ p = ~lwmuf^ (21) 



Q 1 = -^— , 2 = J^i. . (22) 



T -> 73 T , p ->■ 7 f p , 1 -> 73O1 , 2 7 3 2 2 . (23) 



(24) 



At high temperature, the solution is simply the Reissner-Nordstrom AdS solution. The scalar 
profile X(z) being zero means that there is no condensate, i.e., the operators 0\ and O2 vanish. 
The solution is given by taking (we restore dimensionful z for now): 



A(z) = 0, X (z) = 0, * (z) = ^ i_A , /(z) = 1 + g2* * (i + Q 2) . (25) 



1 > /(«) = i + Q -r --3 

So we read off the temperature, chemical potential, and density as: 

T= a (3-Q), n= — , p= -3 - . (26) 
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2.2 Intermediate Temperature Phase 

At low enough temperatures, a new type of solution is available that is a charged black hole with 
a non-zero charged scalar profile. We have non-zero X(z) and x( z ): ancj - we require either Ai = 
or A2 = 0, corresponding to having either O2 or 0\ turned on respectively. 

The solutions can only be exhibited numerically, and we display the resulting plots of 
temperature versus operator vev for each case of O2 and 0\ in figure [2] Below a critical tempera- 
ture T c , this type of solution is thermodynamically favored over the Reissner-Nordstrom case, and 
represents the superconducting phase, with non-zero condensate. These thermodynamics will be 
reviewed in the next section. 



0i 



Vp~(16* G4) 1 ' 




R 1 ' 2 



0.02 0.04 0.06 0.08 0.10 0.12 ^ (16 ttG 4 ) 1/4 
(a) Ox 




ioor r 1 



0.05 0.10 0.15 0.20 0.25 0.30 0.35 
(b) G> 2 



(167rG 4 ) 1/4 



Figure 2: Plots of operator versus temperature for scalar charged black hole solutions with either 0\ or C 2 
non-zero. The vertical dotted line on the O2 plot denotes the transition temperature. See text. 



2.3 Zero Temperature Phase 

The zero temperature solution is an RG flow between two AdS spaces [2T], In the IR, the fields 
have the behaviour: 

\(S) = ln(2 + V5)+\ 1 r a +... , ^)=^~z-P+... , /(*) = |+..., x (z) = Xo+... , (27) 
where 

a= /303_3 /247_1 

01 V 28 2 ' P V 28 2 ' { ' 
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2.4 Thermodynamics 

The on-shell regularized action is given by: 



I = I EH + Id + Ict , with 

Ia = -jr^r [ d 3 xVhK , 

07TG4 J 

Ic - = li/^l-iek/^l^ 2 ' < 29 » 

where we have used the equations of motion to simplify the on-shell action. The quantity V is the 
volume of the M 2 upon which the field theory resides. Putting everything in we get |21j : 

" (/ 3 - 4A 1 A 2 ) > ) = PVG , (30) 



where Q is the Gibbs energy density. So we define the free energy density as: 

T = G + pn = X ^ G R (h - e*°^M>i) . (31) 

The Reissner-Nordstrom free energy density is given by: 

1 



Frn 



16ttG 4 R 

where (zq,Q) are found by solving: 



^(1 + Q 2 )+A), (32) 



-^(3-g 2 )=i?T, 2 ®=-i; 1 e^. (33) 

We define the difference of the free energy densities: 

AT = T RN -T . (34) 

We show the free energy density difference as function of temperature in figure [3| for each case of 0\ 
and 02- When A J 7 > 0, there is a transition from Reissner-Nordstrom to the black hole with scalar 
profile, representing the superconducting phase. This defines the phase transition temperature, T c . 

Note that in the case of 02, it is the upper branch (the choice with higher vev for 02) that 
is favoured. Here, in contrast to the 0\ case where the vev rises continuously from zero at T c , the 
02 operator has a jump in the vev at T c . The physics of these cases is described more in ref.|21|. 

It is worth noting here that while the phase structure in the 02 case seems strikingly different 
from that of 0i, we are aware of another model in the literature that shows how the two can be 
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connected, although using an apparently different mechanism. In refs. [241 125j . the introduction 
of a background current in the M 2 can continuously deform the solution space of the 0\ into that 
of the 02 case. Here, instead of a current, the unusual behaviour of the metric function f(z) is 
responsible for the multivaluedness of available scalar black holes for some ranges of temperatures. 
In either way of thinking about it, there is a new length scale in the theory that manifests itself as 
a finite value in the jump of the free energy and of O2 as the transition temperature is crossed. As 
we will see, the entanglement entropy will be able to detect this new length scale. 
16 ttG 4 rat 




(a) 0i (b) 2 

Figure 3: Free energy density difference. When AJ 7 > 0, the superconductor is thermodynamically favored. 
This occurs at T c « 0.1199^^^ for the O x case and 100T C » 0.3638 ^l^f^ for the 2 case. 



3 Entanglement Entropy 

With a complete supergravity solution in hand, we are ready to study the entanglement entropy at 
each temperature and study the physics.. We use a strip geometry, as outlined in the introduction. 
We choose the following embedding: 

e = x , e = y , z= z(x) . (35) 

The resulting entanglement entropy is given by: 

r t/2 r 2 ( z'(x) 2 \ 1/2 

4G 4 5 = L / dx\ 1 + . (36) 

J-t/2 z V /(*) / 

The extremization problem has a constant of the motion given by: 

111 

(37) 



z 2 z 2 [\ , VW 



where is the location in z of the bottom of the extremal surface. This allows us to write the 
entanglement entropy as: 



4G 4 S = 2LR 2 [ * dz Z \ - 1 = 2LR 2 ( s + 1 ) , (38) 

' Z 2 ^4 - 2 4) f{z) \ • ) 



C 
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where s has dimensions of inverse length with no divergences. The length I is given in terms of z*: 



y/(z*-Z*)f(z) 



(39) 



Under this scaling of equation (16), £ and s scale as: 



t -> 7 3 l£ > « ~» 73^ , 

so we will focus on the following dimensionless quantities: 

s 



VP?, 



(40) 



(41) 



We are now ready to explore the results. 
3.1 0\ Superconductor 

We show in figure [4] the results for s obtained by fixing the temperature and varying £, the width 
of the strip. Larger £ probes more deeply into the infra-red. We show cases with temperature 
below the transition temperature T c . (At the transition temperature, the curves for the Reissner- 
Nortstrom case and the scalar charged black hole (representing the superconducting phase) are 
identical.) In all cases, the curves go linearly with I for large I as is expected from the area law. 
As the temperature is lowered, the slope of the curve for large £ (still linear) is smaller for the 
superconducting background, and continues to flatten out as we approach zero temperature. This 
is expected since for T = 0, the background is an RG flow from one AdS vacuum to another, and 
this flattening out of the entanglement entropy at some finite value was observed in our studies of 
entanglement entropy along RG flow presented in ref. |26j . 
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Figure 4: Entanglement entropy vs. strip width £ for the 0\ case. The solid blue curve is the superconductor 
solution, and the red dashed curve is the Reissner-Nordstrom solution. 
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The fact that for a particular £, the superconducting solution exhibits a lower entropy than 
the Reissner-Nordstrom solution fits with the expectation that degrees of freedom have condensed 
and so there should be fewer of them. 

It is instructive to slice the data differently, fixing a strip width I and studying how the 
entropy of the system changes with temperature. We present this in figure [5] In reading the 
figure, determine the physical curve by always choosing the point of lowest entropy at a given T . A 
discontinuity in the slope of the decreasing entanglement entropy can be observed at the transition 
temperature (indicated by the vertical dotted line), showing its utility as an independent probe of 
the phase structure of the superconductor. 

The slope may be thought of as a sort of response function characterizing the system, 
roughly analogous to a specific heat. It is natural for it to be positive, since increasing temperature 
should indeed promote entanglement entropy. A discontinuous change in the slope at the transition 
temperature T c signifies a significant reorganization of the degrees of freedom of the system. Since 
there is a condensate generated, it is also to be expected that there is a reduction in the number 
of degrees of freedom as well, although precisely at T = T c , the condensate value only just begins 
to rise from zero. We will see something more dramatic in the O2 case, by way of contrast. 

s R 1/2 



(16 ttG 4 ) 



1/4 
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(16 ttG 4 ) 
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Figure 5: The entanglement entropy in the 0\ function of temperature, for fixed £. (We choose 

^fp{l&irGi) 1 / i Rr 1 / 2 l/2 = 2.5) The solid blue curve is from the superconductor solutions, while the red 
dashed curve is from the Reissner-Nordstrom solutions. Trace the physical curve by always choosing the 
lowest entropy at a given T. There is a discontinuity in the slope of the decreasing entanglement entropy 
at the transition temperature T c , indicated by the vertical dotted line. (While we do not plot all the 
superconductor points, due to lack of numerical control at low temperature, we display the zero temperature 
solution, since the solution is known exactly there.) 
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3.2 2 Superconductor 

We show in figure [6] the results for s obtained by fixing the temperature and varying £, the width 
of the strip. As before, larger £ probes more deeply into the infra-red. We again show cases 
with temperature below the transition temperature T c , although in this case, at T c , the Reissner- 
Nordstrom curve lies above that of the superconductor curve. This will mean a discontinuous 
jump in the value of the entropy, as we will see below, in contrast to the 0\ case. For all cases, 
the typical behavior of the curve initially resembles that of the 0\ case, in that the curve of the 
superconducting solution has a lower slope. However, for T < T c , we find multiple solutions for a 
given range of strip widths £, which form a swallowtail shape in our curves. This means there is a 
kink in the physical curve for the entanglement entropy, since we must choose the lowest value. For 
future reference, we will refer to the two parts of such kinked (s, £) curves as the "small £" branch 
and the "large £" branch, respectively, going from small to large £. 

The kink moves to lower strip widths £ as the temperature is decreased (for fixed charge 
density), and persists at zero temperature. Another interesting point is that the leveling off of the 
curve at zero temperature occurs at a positive finite value. This is interesting since in ref. [26J, 
only negative finite values were observed. Furthermore, the 0\ case of the previous subsection also 
exhibits a negative finite value. In ref. [26], we predicted from our sharp domain wall analysis that 
a positive finite value would develop if the domain wall was sufficiently sharp and near the UV. We 
will explore this in subsection |3.4| 




Figure 6: The entanglement entropy for the O2 case. The solid blue curve is the superconductor solution, 
red dashed curve is the Reissner-Nordstrom solution. There is a multivaluedness at finite £ in the shape of 
a swallowtail curve. At lower T, it occurs at lower £, and at T — the swallowtail region persists at finite i. 

It is also interesting to track, as a function of temperature, the strip width value, £k, at 
which the kink appears in the entanglement entropy. We do this in figure [7J (Note also that, for 
fixed temperature, £k oc p^ 1 / 2 .) From this curve we can read off an interesting piece of information. 
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If we pick a specific strip width value £ between an £ m i n ~ 2.27 and an £ m8iX — 2.70 we can read 
off a specific temperature at which the entanglement entropy is crossing over from the "small £" 
branch of an (s, £) curve to the "large £" branch of the (s,£) curve, the branches being separated 
by a kink. The significance of this temperature will be apparent shortly. 
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Figure 7: The position of the kink in the entanglement entropy as a function of temperature. The red 
dashed line is the best fit curve given by ae bx + c, with a = 0.00553, b = 12.9848, and c = 2.25678. The 
curve exists between 3-11 ^min — 2.27 &nd cin ^max 

~ 2.70, discussed in the text. The vertical line indicates 

the position of T c , above which there is no kink. 



A swallowtail multivaluedness in the entanglement entropy, showing multiple extremal sur- 
face solutions at a given £, was first observed and characterized in our studies, presented in ref. |27j . 
of the evolution of the entanglement entropy after a quenching proces^} We will discuss the origins 
of the multivaluedness of the present case in the next subsection. 

Physically, the appearance of a kink in the entanglement entropy as we go to larger £ can 
be attributed to the existence of a new scale in the theory, and the entanglement entropy is a good 
probe of its presence. In contrast to the 0\ case, the transition at T c was associated with a finite 
jump in the free energy J 7 , and also a jump in the the vev of the operator 02- At every subsequent 
temperature T < T c at which this phase is favoured, the system reflects the presence of this scale. 
In terms of the operator 02, there is a length scale that is to be interpreted as a finite correlation 
length, £, inversely related to the vev of 02 at T c . For very small £, the entanglement entropy will 
not be sensitive to £, but when £ becomes comparable to £, correlations between quanta on these 
scales will effectively reduce the number of effective degrees of freedom, reducing the contributions 
to the entanglement entropy, resulting in a kink to change of the slope of the (s, £) curve for larger £. 

As was done for the case of 0i, it is instructive to study the entropy for fixed £ as the 
temperature varies. There are in fact three distinct situations, giving rise to three different types 
of curve. The distinguishing issue is whether the choice of fixed £ can ever become a kink value, 



Refs. |28l 129] have also since observed this phenomenon. 
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at some temperature. As we saw, figure [7] has the answer to this. 

The first case is to have a fixed £ that is greater than ^ max - Then, at successively lower 
temperatures than T c , the entropy will always come from points on "large C branches of (s, £) 
curves, since the kink moves to smaller £ as T is reduced. We show the resulting (s, T) type of 
curve in figure [8j In this and the next two figures, the red curve is always favoured for T > T c , 
indicated by the vertical dotted line. Below T c , one should determine the physical curve by always 
choosing the point of lowest entropy at a given T. Notice how the region of negative slope is nicely 
avoided. This will be the case in our subsequent curves as well. 

s R 1 ' 2 
JJ (16 ttG 4 ) 1/4 




Figure 8: The entanglement entropy in the O2 case, as a function of temperature, for fixed I. (We choose 
^(167^4) 1 l i R- 1 l 2 l/2 = 3 here.) The solid blue curve is from the superconductor solutions, red dashed 
curve (decreasing in slope, but only slowly on this scale) is from the Reissner-Nordstrom solutions. Trace 
the physical curve by choosing the red curve for T > T c , indicated by the vertical dotted line, and for T < T c , 
always choosing the lowest entropy. There is therefore a discontinuous jump in s and its slope at T c . (While 
we do not plot all the superconductor points, due to lack of numerical control at low temperature, we display 
the zero temperature solution, since the solution is known exactly there.) 



The second case is to have a fixed £ that is smaller than £ m i n . There is no temperature at 
which this would become a kink value, since although the kink moves to smaller £ with smaller T, 
it stops at £ m i n at T = 0. So the contributions to the entropy come entirely from points on "small 
£" branches of (s,£) curves. We show the resulting (s,T) type of curve in figure [9j 

The final case is to have a fixed £ such that £ m \ n < £ < £ ma ,x- Then, as we reduce the 
temperature from T c , entropy contributions are from "large £" branches until a temperature is 
reached such that our chosen £ is a kink value £k- This temperature can be read off from figure [7| 
For lower temperatures, the entropy will be from points on "small £" branches. Consequently, the 
superconductor phase of the (s, T) curve in this case will be a combination of two types of curve, 
connected by a new discontinuity in the derivative where they join. We show the resulting (s, T) 
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Figure 9: The entanglement entropy in the 02 function of temperature, for fixed £. (We choose 

v /p(167rG 4 ) 1/4 i?" 1/2 ^/2 = 2.2 here.) The solid blue curve is from the superconductor solutions, red dashed 
curve (decreasing in slope, but only slowly on this scale) is from the Reissner-Nordstrom solutions. Trace 
the physical curve by choosing the red curve for T > T c , indicated by the vertical dotted line, and for T < T c , 
always choosing the lowest entropy. There is therefore a discontinuous jump in s and its slope at T c . (While 
we do not plot all the superconductor points, due to lack of numerical control at low temperature, we display 
the zero temperature solution, since the solution is known exactly there.) 



type of curve in figure 10 



In all the curves in figures [8j [9] and figure [TOj we see that in addition to the slope having a 
discontinuity at the transition temperature T c (shown by the vertical dotted line) , the value of the 
entropy drops discontinuously as well, cis we earlier anticipated. Then, when -^min 
have the additional feature of a discontinuity in the slope at some lower temperature, generated (as 
discussed above) by the sensitivity to the length scale £. Note that for fixed I with a value close 
to the edges (but outside) of this interval, the resulting (s, T) curves will appear to have a locally 
smoothed out discontinuity. The cases displayed here are far away enough from the interval that 
the smoothing is spread out. 

3.3 Multivaluedness of the Entanglement Entropy 

As stated already, we've seen a multivaluedness of the entanglement entropy before, in a study of its 
evolution after a quench |27j . This resulted in a kink representing the change in the saturation rate 
of the entropy as it evolved. To understand the reason for the appearance of the multivaluedness in 



this case, we study the behavior of the function f(z), which we show examples of in figure 12 (we 



show the behavior for the 0\ function in figure 11 for comparison purposes). Slightly before the 



rightmost point in the vev curve of figure 2(b)[ f(z) develops two new extrema, a minimum and a 



maximum. The maximum grows higher in value and the minimum becomes sharper (i.e. greater 
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^7 (16 7rG 4 ) 1/4 
0.55 




100 T R l/2 



0.5 JJ (l6nG 4 ) 1/4 



Figure 10: The entanglement entropy in the C 2 function of temperature, for fixed £. (We choose 

y /p(16TrG 4 ) 1 ^R- 1 ^£/2 = 2.5 here.) The solid blue curve is from the superconductor solutions, red dashed 
curve (decreasing in slope, but only slowly on this scale) is from the Reissner-Nordstrom solutions. Trace 
the physical curve by choosing the red curve for T > T c , indicated by the vertical dotted line, and for T < T c , 
always choosing the lowest entropy. There is therefore a discontinuous jump in s and its slope at T c . There 
is an additional jump in the slope at a lower temperature. (While we do not plot all the superconductor 
points, due to lack of numerical control at low temperature, we display the zero temperature solution, since 
the solution is known exactly there.) 



second derivative) as the temperature decreases. Both extrema move toward the AdS boundary at 
z = (the UV) as the temperature decreases. 



We show the behavior of the function f(z) at zero temperature in figure 12(b) As for 
the finite (but low) temperature case, the function has a minimum. This minimum is now a finite 
distance from the UV boundary, and the maximum we saw at finite temperature has now smoothed 
out to a constant as we go toward the IR. It is the non-monotonic behavior of f(z) that generates 



the multivaluedness, as illustrated in figure 13 When z*, the location of the lowest point of the 
extremal surface, is in the neighborhood of the minimum of f(z), the entanglement entropy becomes 
multivalued. Only when relatively far from the minimum does the entanglement entropy become 
single-valued again. 

3.4 Domain Wall Behaviour 

The non-monotonic behavior of f(z) suggests that the behavior of the domain wall that interpolates 
between the two AdS vacua (at zero temperature) is also not monotonic. In addition, we'd also 
like to compare the domain wall features of the cases 0\ and 02, to see if we can understand the 



differences between the large i saturation values of the entanglement entropy observed in sections 3.1 



and 3.2 (negative versus positive) in the terms discussed in our RG flow studies of ref. |26j. To 
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Figure 11: The metric function f(z) for the 0\ superconductor. For (a), red dashed ( (y^g ) 1/2 ~jp ~ 

0.02990), solid blue 



0.11993), green dot-dashed ^ u/2 4= = 0.10309), cyan dotted (jt-^ 
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Figure 12: The metric function f{z) for the C 2 superconductor. For (a), red dashed ( — 
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study this, we first write the metric in the following form: 

ds 2 = _ e 2A{r) e - X (z) +X (0) dt 2 + + ^2 (42) 

z(r)^ 

The function ^4(r) encodes the domain wall, and the coordinate r can be determined via: 

rW - r °=-£iW)- (43) 

Note that technically ro is actually at — oo since that corresponds to the IR of the theory. So to 
circumvent this issue, we can define a variable f such that: 

and simply shift our result for A(r) by a constant Aq such that the ratio (A(f) — Ao)/r does not 



diverge at r = 0. We show the behavior of A{f) in figure 14 Indeed, as expected, it is non- 
monotonic for the 02 case, and indeed the domain wall is much sharper for this case than it is for 
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(a) f(z) (b) Entanglement Entropy 



Figure 13: A comparison of where the minimal surface corresponding to the entanglement entropy ends 
and the multivaluedness of the entanglement entropy. The colored and numbered dots correspond to where 
the surface used to calculated the entanglement entropy ends in the bulk. 



the 0\ case, confirming our observations made in section 3.2 and ref. [26] about the saturation of 
the entanglement entropy at large i at T = 0. The multivaluedness feeds nicely into the swallowtail 
structure, as we saw in the previous section by looking directly at f(z). 



A(f)-A 



A(r)-A 



4 6 

(a) 




(b) 



Figure 14: The domain wall for (a) the 0\ case and (b) the 02 case. The UV is at f oo, and the results 
have been shifted such that the IR results are at finite values of the radial coordinate. 



4 Concluding Remarks 

Since we have carefully unpacked and discussed our results during our presentation of them, we 
will be brief in this section. 

We have presented a study, using holography, of the entanglement entropy of a certain type 
of strongly coupled superconductor. Since the background is fully backreacted and highly stable 
(in the sense outlined in footnote 3), we can be confident that the results are robust. 

This is the first such study of its type, and the results may well be of interest beyond the 
confines of holography, since it is of interest in the condensed matter physics community to use 
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entanglement entropy as a probe of new physics of experimental relevance. Indeed, we have found 
that the entanglement entropy is a very sharp probe of the physics at the transition temperature T c , 
the ground state of the system at T = 0, and also at intermediate temperatures, where, in mapping 
out the full temperature range, we identified various kinds of novel transition in the entropy. Some 
of the novel physics (arising from multiple extremal solutions for the entropy at a given point 
in parameter space) recalled phenomena observed in our earlier studies of the time evolution of 
entanglement entropy [27] . 

The study presented here also served as another holographic example of the behaviour of 
entanglement entropy along an RG flow, which we studied in ref. [26], and in fact we were able to 
confirm some more of our predictions from that paper using phenomena observed here. 

We expect that this is just the beginning of a series of fruitful investigations of this type, 
shedding more light on a wide variety of strongly coupled quantum phenomena of interest. 
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